Lifshitz Hydrodynamics 
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We construct the hydrodynamics of quantum field theories with a Lifshitz scaling symmetry. New 
transport coefficients are allowed by the absence of boost invariance, however, only one is compatible 
with a local increase of the entropy density. The formulation is applicable, in general, to fluids with 
an explicit breaking of boost symmetry. We use a Drude model of a strange metal to study the 
physical effects of the new transport coefficient. It can be measured using electric fields with non- 
zero gradients, or via the heat production when an external force is turned on. Scaling arguments 
fix the resistivity to be linear in the temperature. 



Introduction.- Quantum critical points [lj, [2| are 
believed to be responsible for the strange metal be- 
haviour observed in a variety of materials including 
heavy-fermion compounds and high T c superconductors. 
Although the quantum critical point is strictly defined at 
zero temperature, its properties affect a large region of 
the phase diagram. 

As in ordinary critical points [3j , there is a scaling sym- 
metry at the critical point acting on the space-time co- 
ordinates as 



t ->• fl z t, 



c' — > fla 



= I,.. A 



(1) 



We assume that at the relevant length scales the system 
is isotropic and homogeneous, thus rotation and trans- 
lations (both in space and time) are symmetries. All 
together this is the 'Lifshitz symmetry group', inspired 
by the scaling in Lifshitz points [J, [5[ . 

We use symmetries and the laws of thermodynamics 
as guiding principles in order to write the hydrodynamic 
equations of quantum critical points at finite tempera- 
ture. We expect to have a universal expansion in deriva- 
tives of the velocity and temperature with some coef- 
ficients, whose precise values depend on the microscopic 
theory. Fluid dynamics have been shown in many cases to 
describe strongly correlated critical systems, such as con- 
formal field theories at finite_temperature |6[ , fermions at 
unitarity [7| and graphene [8l4l0l|. 

An important aspect of theories with Lifshitz scaling 
symmetry is the absence of boost invariance. The energy- 
momentum tensor, which is symmetric for theories with 
Lorentz invariance, can receive now asymmetric contri- 
butions, and in principle many new transport coefficients 
can appear. We use the second law of thermodynamics 
in its local form in order to constrain the transport co- 
efficients and find that there is a unique one. We then 
take the non-relativistic limit and identify the dissipative 
terms related to the breaking of boost invariance. The 
results are applicable in general to fluids with an explicit 
breaking of the boost symmetry, not necessarily with Lif- 
shitz scaling. 

Since Lorentz or Galilean invariance are fundamental 
symmetries in any microscopic theory, for the case of 



a strange metal we interpret its absence in the fluid as 
the effect of an external medium through which the fluid 
moves. Following this philosophy we study a simple gen- 
eralization of Drude's model, where rather than a weakly 
interacting gas the collective behaviour of the electrons 
is an almost perfect fluid. The effect of the medium is a 
drag term and the new dissipative term. A similar model 
for graphene was considered in ll[ , but to our knowledge 



the effects due to the breaking of boost invariance have 
never been considered before. We propose different mea- 
surements of the conductivity and the heat production 
that can be used to determine the new transport coeffi- 
cient. 

Interestingly, we also find that scaling arguments fix 
the resistivity to depend linearly on the temperature 
p xx ~ T, which is the observed behaviour in experi- 
ments 12], under the reasonable assumption that the 
dependence on the density is linear. The result seems 
to be universal, in the sense that it is independent of 
the number of dimensions or the value of the dynamical 
exponent. 

Hydrodynamics.- The generators of Lifshitz sym- 
metry are time translation Pq — dt, spatial translations 
Pi = di, the scaling transformation D = —ztdt — x l di 
and rotations. The subalgebra involving D, Pi and Pq 
has commutation relations 



[D, P % ] = P % , [D, P ] = zP Q 



(2) 



In a field theory the scaling symmetry is manifested as 
a Ward identity involving the components of the energy- 
momentum tensor 



zT u , 



S\T) = . 



(3) 



At finite temperature T° = —e, T 1 ^ = p5*j, leading to 
the equation of state 



ze = dp. 



(4) 



The Lifshitz algebra can be generalized for constant ve- 
locities u M , u^u^ — r\ ilv u^u 1> = — 1. We define the gener- 
ators 



PU = u"0„, P£ = P;d v , D = zx»u^P n -x fl P t 



(5) 



Where P ' " = 8" + u^u v . Then, the momentum opera- 
tors commute among themselves and 



[D,P"]=zP\ [D,P^] = P^. 
The Ward identity associated to D becomes 



zTt/, 



TtP„ 







(6) 



(7) 



It coincides with (|3|) only when z = 1, but leads to the 
equation of state (2)) for any velocity. If one insists on 
preserving the Ward identity in its original form ([3]) , then 
the resulting equation of state will be velocity dependent 



ze - 



dp+(z~ l)(e + p)uj 







(8) 



We do not study this possibility, rather we take (0) as 
the correct form after a change of coordinates to a frame 
where the fluid is moving. Our choice is supported by the- 
analysis of the model proposed in [13[ as a gravitational 
dual to Lifshitz points, and we will detail the calculation 
elsewhere. 

The conservation of the energy-momentum tensor 
determines the hydrodynamic equations d^T^ = 0. 
Lorentz symmetry forces the energy-momentum tensor 
to be symmetric, but in other cases, as in Lifshitz the- 
ories, this condition can be relaxed. This allows many 
new terms in the hydrodynamic energy-momentum ten- 
sor, but as usual there are ambiguities in the definition 
of the velocity and the constitutive relations. In order to 
fix them, we impose the Landau frame condition 



r p^ l 



-ew 



(9) 



Then, the generalized form of the energy- momentum ten- 
sor is 



T^ ={e + p)vPu v +prT 



Av) 



M v 



[va] 



Jmo-] 



+ nY i + {u^T i +^^P) 



(10) 



The first line is the ideal part of the energy-momentum 
tensor, while in the second line we write the dissipative 
corrections. 7r^ is antisymmetric in its two indices and 
its is symmetric. In addition, the Landau frame condi- 
tion imposes ij m„ = 0. In a theory with rotational 

[ill 

invariancc n^ = 0. There may be other new transport 
coefficients in a theory with more conserved charges, how- 
ever for the example of a single conserved global current 
we do not find any to first order in derivatives. 

Assuming Lifshitz symmetry, the scaling dimensions of 
temperature, energy, pressure and velocities are: 



[T] 



[p} = z + d, 







(11) 



Entropy current.- The new terms should be compat- 
ible with the laws of thermodynamics, in particular with 
the second law. Its local form in terms of the divergence 
of the entropy current is d^j£ > 0. The divergence of 



the entropy current can be derived from the conserva- 
tion equation 



= d^Uu = -Tdr(8U») 



(12) 



The dots denote contributions originating in symmetric 
terms in the energy-momentum tensor, to first order in 
derivative corrections they will simply be the shear and 
bulk viscosity contributions, which are manifestly posi- 
tive for positive values of the transport coefficients. In 
order for the antisymmetric contribution to be positive, 
we should be able to write it as a sum of squares. This 
is possible only if 



itAlna] = at(na)(di t j.Uv] - u[ p ,u a d a u <T ]) 



(13) 



where a{jiv) = a^is/j,), a.(p,p) = are new transport 
coefficients. The second law demands 



a(0i) > 0, a(ij) < 



(14) 



The coefficients a(0i) break boost invariancc and a(ij) 
rotational invariance. If only boost invariance is broken, 
then a(0i) = a, Vz. For a theory with Lifshitz symmetry 
the scaling dimension is [a] = d, which determines the 
temperature dependence of the new transport coefficient 
to be 



T* 



(15) 



Kubo formula for isotropic fluids.- Under the trans- 
formation x 1 — > x l + £ l , the expectation value of the 
energy-momentum tensor is shifted, to leading order, by 

(T^(x))= [ d d+1 x'(T^{x)T a i{x'))d a ^{x') . (16) 



doC can be identified with a 



The time derivative v % 
velocity. Then, 



n m = </ T [o^ = f d d+i x , ^ T [oi] T o^ v . (17) 



We have used that (T^T^) = by rotational invari- 
ance. Let us Fourier transform to frequency w and mo- 
mentum q space the velocity and the correlation func- 
tions. Then, the first order transport coefficient is given 
by the Kubo formula 



Q 



-liml% 



Im/T^T Q 3\(uj,q = 0) . (18) 



Non-relativistic limit.- The relativistic velocity is 
w/' = (l,P l )/^l — /3 2 , where f? — v l /c. In the non- 
relativistic limit Pi — > 0, the pressure is not affected while 
the relativistic energy is expanded in terms of the mass 
density p and the internal energy U as 



2 

c p 



2 



U 



(19) 



The relativistic hydrodynamic equations reduce to the 
non-relativistic form 



d t p + d i (pv i )=0, (20) 

d t U + di (Uv l ) + pdy 

= ~v ij cn J + ^(diV i ) 2 + ~(Vi) 2 , (21) 

dtOw*) + d j (pv i v i ) + 8 i p 

= dj (w ij + ^S ij d k v k 

+ d t (aVX) + d j (| (y*V\ + v'Vi)) . (22) 



The shear tensor is o; 



diVj + djVi - (2/d)Sijd k v k 



While taking the limit, we have absorbed factors of 1/c 
in the shear and bulk viscosities r\ and C and a factor 1/c 2 



in a. The vector V\ is defined as 



VX = 2D t v l + w v v i , (23) 

where D t = dt + v l di is the material derivative and uiij — 
2d\iVj] is the tensor dual to the vorticity. The first term in 
V\ is proportional to the relative acceleration of the fluid, 
while the second term is proportional to an acceleration 
due to the Coriolis effect. Similarly to the viscosities, the 
coefficient a determines the dissipation that is produced 
in the fluid when the motion is not inertial. 

The non-relativistic limit with p ^ does not allow 
Lifshitz (or for that matter conformal) scaling in the 
original relativistic theory. However, one can have a non- 
relativistic version of Lifshitz scaling. Under a space-time 
diffeomorphism 

*-►* + £*, x i ^x i + C , (24) 

the partition function of the theory will change as 

6logZ= [dtdtxi-d^jg + dnCTl) . (25) 

The Lifshitz equation of state is recovered if the theory 
has a symmetry 



£* = zt, f = x i + - — 2 -v l t 



(26) 



This is a combination of a scaling transformation (JTJ) and 
a change of frame. The Ward identity becomes 







-"-.; 



t+J2 T '\ 



z-2 



v l T\ = -zU + dp . (27) 



Note, that in the relativistic case the analogous term in 
(|4|) that modifies the equation of state is proportional to 
z — 1 rather than z — 2. The reason is that, for 2 = 1 
in the relativistic case and z — 2 in the non-relativistic, 
the Lifshitz symmetry group can be extended to the con- 
formal groups that include boost transformations. The 
scaling Ward identity is then determined by a transfor- 
mation that is independent of the velocity. 

Conductivity.- Consider a charged relativistic fluid 
moving through a static medium with an external elec- 
tromagnetic field applied to it. The hydrodynamic equa- 
tions are 

<V = 0, d^ v = F^Jf, - XcPfJv , (29) 



where P^ v 



rf v + tH u ,t li = (1,0), is the projector on 
the directions transverse to static medium. The coeffi- 
cient A determines the drag force that the fluid feels as 
it moves through the medium. 

Let us define E, = Fqi and consider the case where 
the magnetic field is zero and d^Ei = 0. Projecting the 
equations of the energy-momentum tensor with t M and 



P t , we find 



d^ = J l Ei 



dpTi* = J°E l - XcJ 1 



(30) 



We are interested in describing a steady state where the 
fluid has been accelerated by the electric field, increas- 
ing the current until the drag force is large enough to 
compensate for it. We assume that the flow does not 
change, but some scalar quantities like the energy can 
change with time. In order to simplify the calculation 
we will take the non-relativistic limit and consider only 
an incompressible fluid diV 1 = 0. The fluid motion is 
described by the Navier-Stokes equations 

pv k d k v i + d l p (31) 

= pE l - Xpv 1 + ?7VV + ^dj {{v 3 a lk + uV fc ) v k ) . 

We can solve this equation order by order in derivatives, 
keeping the pressure constant d l p = 0. To leading order 
the current satisfies Ohm's law 



J 1 = pv l ~?-E l 
A 



(32) 



and the conductivity is simply er^ = p/X6ij. At higher 
orders in derivatives we find the following corrections for 
a divcrgcnceless electric field E x (y) 



,{E X ) = £ 



1 



^ (■*** + &<«* 



p\E x 



(33) 



As for the relativistic case, imposing a scaling symmetry 
without the change of frame would modify the equation 
of state 



zU + dp 



z~2 



-pv 



(28) 



The conductivity depends on the electric field and its 
gradients. In the case where the electric field is linear 
E x = Eoy/L, the conductivity is simplified to 



1 



pL 2 X 3 



(34) 



Note, that the contribution from the shear viscosity 
dropped. This gives a way to identify the new transport 
coefficient a, it can be measured experimentally as an 
enhancement of the conductivity with the electric field. 
Another simple case is when the electric field takes the 
form E x = Eq cos(y / L) . The contribution of a to the 
conductivity is y dependent 



0xx{y) = 



l- 



v 



aE 2 



ZaEl 



XpL 2 X 3 pL 2 2X 3 pL 2 



(35) 



If we average on the y direction, we find again an en- 
hancement of the conductivity with the electric field 



1 



V 



aE 2 



XpL 2 AX 3 pL 2 



(36) 



Lifshitz scaling.- In a fluid with Lifshitz symmetry 
the scaling dimensions of the hydrodynamic variables are 

[v*] =z-l, [p] = [U] = z + d, [p] = d + 2 - z , (37) 

while the temperature has scaling dimension [T] = z. 
We can determine the scaling dimensions of the trans- 
port coefficients by imposing that all the terms in the 
hydrodynamic equations have the same scaling. We find 



W = z, lv) = [C] = d, [a] = d-2(z-l) 



(38) 



In contrast with a relativistic fluid, the density is ap- 
proximately independent of the temperature. This in- 
troduces an additional scale, and in general the trans- 
port coefficients can be non-trivial functions of the ratio 
t = T i J p. The conductivity will have the following 
temperature dependence 



T- 



-<t(t) 



T 



(39) 



where we assumed a linear dependence on the density as 
obtained from the calculation with the drag term. Note, 
that this predicts a resistivity linear in the temperature 
and independent of the dynamical exponent and the num- 
ber of dimensions (12l |. However, a prediction for the Hall 
angle that matches with the experiment cot 8h ~ T 2 
14j | is not as straightforward. The naive scaling from 
this model would be cot Oh ~ T. The scaling could be 
different if there is a strong temperature dependence of 
the permeability, while the permittivity is approximately 
constant. 

Dissipative effects.- We can introduce constant ho- 
mogeneous forces by electric fields or temperature gradi- 
ents. They will induce an acceleration 



-d'p/p + E l = {s/ P )d l T + E l 



(40) 



We assume dta 1 = 0, djCi 1 = 0. The Navier-Stokes equa- 
tions for homogeneous configurations takes the form 



d t v l - 2(a/p)d$v i + Xv l = a 



(41) 



If we assume that the forces are suddenly switched on at 
t = 0, the evolution of the velocity is determined by this 
equation with the initial conditions v*(t = 0) = 0, 



*><«-«>-£ ^ + i-i 



(42) 



This choice is based on the physical requirement that at 
large times the velocity stays constant. When a — > it 
simply becomes dtv l (t = 0) = a' 1 . 

The heat production rate induced by the force is 



d t U = Xpv 2 + 2a(<9 t v) 2 



(43) 



At late times the system evolves to a steady state con- 
figuration with constant velocity, so the heat production 
rate becomes constant v l = a 1 /X. Subtracting this con- 
tribution for all times, the total heat produced is 



AQ 




(44) 
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